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We derive on the level of quantum optics expressions for the uncertainty of the photocount in a 
multi-mode multi-frequency setup. The result depends on the quantum correlations of the individ- 
ual modes and the frequency spectrum of the radiation, the latter leading to a frequency beating 
sometimes referred to as dynamic laser speckle. When the mode structure of the radiation field is 
disturbed between source and detector, another contribution to the photocount uncertainty referred 
to as static speckle appears. To predict the size of this effect, we present a suitable definition of the 
etendue (or phase space volume) that links the number of modes of a radiation field to macroscopic 
quantities. 



I. INTRODUCTION 

Given some radiation field, a prediction can be made 
about the outcome of an experiment where a suitable 
detector is placed at some position R and probes the 
radiation field for some time T. Even when the radiation 
field is known as precisely as possible, there is still a finite 
uncertainty in the outcome. This statement can also be 
formulated in the inverse way: it is not possible to build 
an apparatus that provides a "better" (=more defined) 
illumination at the point R than this limit. 

We will demonstrate in this paper that the uncertainty 
of the photocount, quantified by its variance, depends on 
both the number of modes, including the energy distri- 
bution among the modes, and on the spectral properties 
and quantum correlations of the radiation field. Another 
contribution to the uncertainty will appear if the radi- 
ation field is known at the source but the mode struc- 
ture is then perturbed in some uncontrolled way before 
it reaches the location of the detector. The latter effect is 
frequently referred to as static speckle 1] named after the 
pattern it creates on a camera. The opposing term "dy- 
namic speckle" unfortunately is used to denote two com- 
pletely different concepts in the literature. Frequently it 
refers to a static speckle pattern that is changing over 
time because the source of the perturbation is moving 3 
but in this paper we refer to "dynamic laser speckle" 3| 
- an effect that is due to the quantum dynamics of the 
radiation field. 

Quantum effects are most prominent on small length 
scales. Modern optical lithography operates precisely in 
this regime, printing structures smaller than the wave- 
length of the light with a precision of a few nanome- 
tres Q . To this end, photo resists have been developed 
that are more sensitive to changes of the light intensity 
than most technical sensors, thereby acting as efficient 
(albeit unintentional) detectors for uncertainties of the 
photocount. Later in this paper we will show that the 
effects of multi-mode multi-frequency radiation are most 
pronounced neither at the limit of large or short times 
but in the intermediate regime. Recent studies have con- 
firmed that this is precisely the regime used in modern 
optical lithography |5|]. 



II. OVERVIEW 

Classically a radiation field is described by its electric 
field E(f, t) as function of position r and time t. In a 
quantum treatment, the electric field E is replaced by a 
suitable operator. In both treatments, the electric field is 
not a perfectly determined quantity, and the maximum 
knowledge possible is contained in the mutual-coherence 
function [6( 

r(rx,r ? 2,ii,i 2 ) = (E*(r 1 ,t 1 )E(f f 2 ,t 2 )) . (1) 

We restrict ourselves to stationary fields such that T does 
not depend on t\ and ti but only on the time-difference 
t\ — 1 2 . Furthermore, for ease of writing we will only treat 
one component E of the electric field E but the extension 
to cover the other components is straight-forward. 

The Fourier transform P(ri,r2,a->) of r(fi,f*2,ti — t%) 
with respect to the time difference is a nonnegative Her- 
mitian operator, and thus possesses an eigenrepresenta- 
tion @| 

f(fi,f 2) w) = a n (u) (ft* (f^ , (J) <)> n (f 2 , w) , a n (u) > . 

n 

(2) 

The eigenvectors (f> n (f,ui) are called the modes of the 
electric field and form a complete orthonormal set. We 
assume that the frequency spectrum is small enough such 
that (j> n (r, ui) is independent of oj. Whenever the fre- 
quency spectrum is equal to the natural linewidth of the 
light source, i.e., the spectrum is due to the finite life- 
time of some excited light-emitting medium such as in a 
laser, this assumption is always fulfilled. For other light 
sources, this assumption can sometimes be problematic 
for cavity-like systems but for open system or in a waveg- 
uide geometry, this is less of an issue. If the condition 
that <ft n (r, w) is independent of ui, should be violated, one 
can split the mode into several discrete modes (one for 
each frequency interval in which the shape of the mode 
can assumed to be constant) and the remainder of this 
paper be applied nonetheless. The electric field can then 
be written as 

E(f,t) = J2"n(t)cl>n(r f ). (3) 

n 



2 



Since <p n is independent of i, the <p n are the modes instead 
of being some arbitrary base of the electric field. Compar- 
ison of Eqs. (|2|) and ((3|) shows that the magnitude of a n 
can be computed from a n but not its phase. A semiclas- 
sical treatment of the photocount statistics would assume 
that a n is fluctuating in time whereas within quantum 
optics, which we will apply in this paper, fluctuations are 
inherent to the operator description. 

This paper is organised as follows. In Sec. IHII we com- 
pute the variance of the photocount on a quantum-optical 
level when the radiation field is completely known. We 
will find that there is a shot noise term, a quantum cor- 
relation term and a term describing the beating of differ- 
ent frequencies. Frequently, by design or unintentionally, 
the mode structure emitted by a known light source is 
completely changed and unknown when it reaches the 
point of the detector. As will be demonstrated in Sec. El 
random-matrix theory allows a compact and exact treat- 
ment of this problem. The resulting "static speckle" be- 
comes the smaller the more modes of the radiation field 
are excited, and the more uniform the energy distribu- 
tion among the modes is. Depending on the "size" of the 
radiation field, there is thus a minimum amount of static 
speckle, which is computed in Sec. IIVI In Sec. IVII we ex- 
tend this question to computing the most likely amount 
of static speckle. Since we will demonstrate this quantity 
to be self-averaging, the computed average is more than 
just an average in that it is characteristic for almost all 
individual speckle values. 



III. QUANTUM THEORY OF 
PHOTODETECTION 



While in "general" electrodynamics the electric field 
E(f, t) can be probed directly, this is not possible in the 
realm of optics as E(f, t) changes too quickly in time 
and space to allow a direct measurement. Rather, the 
radiation field is probed by means of photodetection, i. e., 
by absorbing photons inside some device and counting 
the number of photons absorbed. This principle applies 
to technical machines as well as to the human eye or to 
photo resists. 

On small length and time scales, quantum effects be- 
come important. For this reason, and because all effects 
can then be treated in a more compact way, we will use 
the quantum theory of photodetection in the following. 
This theory [7J, [8( is usually formulated for single-mode 
detection. Since multi-mode fields Q are at the core of 
this paper, we will present an extension of the theory to 
multi-mode photodetection here. We will treat the multi- 
frequency aspect explicitly and not express the different 
frequencies by different modes as is frequently done in 
textbooks. 

We label the modes of the electromagnetic field by 
4> n (r). The annihilation operator associated with this 
mode is a n . Using this notation, the quantum operator 



for the electric field at some position R is given by [6[ 

F(t)=J2<t>n(R)an(t) ■ (4) 

n 

Later it will prove helpful to switch from the time repre- 
sentation a n (t) to the spectral representation a n {ui), 

1 



a n {t) 



/2vr 



a n {u)e~ lut dw 



(5) 



since in time representation, operators taken at different 
times do not necessarily commute whereas they do in the 
spectral representation, 

[a n (oj), <4i(w')] = S nm S(uj - uj') . (6) 

When we label the number of photons in the n-the mode 
with /„, which basically amounts to the total intensity 
of that mode, this gives the expectation value 

(aj,(w)a m (u;)) = 2ir5 nm 8{uj - J)I n G n (w) , (7) 

where G n (ui) is the (normalised) spectrum of the radi- 
ation in the n-th mode, and the prefactor 2w has been 
introduced for later convenience. 

Photodetection within some time interval T is then 
described by the quantity @ 

W = r] [ F^(t)F{t)dt , (8) 
Jo 

where 77 marks the detection efficiency and includes in- 
formation on the size of the detector. This amounts to 
a perturbative description of the interaction between the 
detector and the radiation field. The perturbative ap- 
proach neglects the effect that every detected photon de- 
creases the number of photons remaining in the radiation 
field @ but this is mainly an issue for microcavities where 
only a small number of photons are excited at one point 
in time and the dynamics are then studied. For the pur- 
pose of this paper, this is no relevant restriction. 

The factorial moment := (n(n — 1) • • • (n — k + 1)) 
of the photodetection count distribution is given by 

where the colons denote normal-ordering of the operator 
within, and the brackets (. . .) denote the average which 
has to be taken over both the quantum fluctuations of 
the operators a n and the frequency distribution G n (uj). 
Mean and variance of the photo count then follow from 

n = n^ , varn = ? i« + - [n«] 2 . (10) 
For the first moment, this yields the result 



dtJ2K(R)^m(R)ai(t)a m (t)) = 

nm 

± dt ff d Wl , 2 ^^(i?)0 m (i?)(at( Wl )a m (u; 2 )>e i (- 1 - 

nm 

= v [ dt /d^V|0„(i?)i 2 /„G„M = Tr7y>„(i?)i 2 /, 

Jo J 



u 2 )t 



(11) 
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This is the same result as expected by classical theory 
since \(f> n (R)\ 2 I n corresponds to the intensity of the n-th 

I 



mode at position R. 

The second factorial moment is 



n^=r, 2 ( 



(1*1,2 (^)^»2 {RWn 3 {R)<t> ni (-R)a^ {h)a n2 (h)al 3 (t 2 )a n4 (t 2 ) 



a i (uj 1 -u 2 )t 1 + (uj 3 -uj 4 )t2 



(12) 



Taking the average over the quantum mechanical expectation operators gives nonzero contributions for three cases of 
the indices n: ri\ = n 2 = = 714, n\ = n 2 7^ ri3 = «4, and ri\ = 714 7^ n 2 = n^, thus 



,(2) 



Ho dtl 2 1 dUJl ' 3 ^ l^i^l^^a^l^^CwOOni^Oo^Ma^Cws)) 



4tt 2 



4tt 2 



dwi,s 



di 



1,2 £ l^)! 2 !^)! 2 ^^)^^ . (13) 



r 



This equation can be simplified by inserting the expecta- 
tion values from Eq. (JT]) . The double integral over ti and 
t 2 in the last line of Eq. (1 13[) can be reduced to a single 
integral via 



f dt li2 f(h-t 2 )= f dtf(t)[T-\t\ 

J-T 



(14) 



The remaining two integrations over lo\ and to 2 in that 
line each yield the Fourier transform G(t) of G(cu), 



G(t) 
This thus gives 



1 



(15) 



,(2) = 



ifT 2 
1^2" 



ifT 2 £ \<t> ni (R)\ 2 \<l>n 3 (R)\ 2 InJn 
2 i-T 



J dt H \KAR)\ 2 \^{R)\ 2 

xI ni I n2 G ni (t)G*(t)[T-\t\_ 



In J _t 



(16) 



The summation over n\ 7^ 77,3 is rather inconvenient. 
Thus, we include the missing terms ri\ = in that sum- 
mation, yielding the value [n*- 1 -*] 2 , and correct for this by 



subtracting in the first summation the terms just added. 
There now appears the quantity 

full := -±3 J M\al HlVM] 2 ) - (4 Ma„H) 2 , 

that quantifies the deviation of the radiation in the n- 
th mode from a coherent state. For coherent radiation, 
/„ = since then ([a* (w)] 2 [a„(w)] 2 ) = {al(uj)a n (uj)) 2 
whereas for any Gaussian light source, in particular any 
object emitting thermal radiation, ([aj l (w)] 2 [a„(u;)] 2 ) = 
2(a^(w)a n (w)) 2 , and f n = 1. For classical radiation, 
fn > whereas for certain nonclassical radiation /„ < 
is possible. 

Collecting results, using Eq. (fTO)) . gives the variance of 
the photocount 



varn = n + rfT 2 V|0 n (i?)| 4 / 2 /„+ 



rf_ 
2tt 



£ \cb ni (R)\ 2 \4>n 2 (R)\ 2 I ni In 2 

x / T diG ni (t)G; 2 W[T-|i|]. (18) 



Equation (fTg|) constitutes a core result of this paper. 
It gives the general expression for the noise that is intrin- 
sic to any measurement of a multi-mode multi-frequency 
radiation field over some time T. The first term is the 
well-known shot noise term that appears naturally in the 
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quantum treatment: electric field and photocount are 
quantum-optically related via the factorial moment, see 
Eq. ((9J, whereas in classical optics the relation is via the 
"plain" moment, with the difference between factorial 
and plain moments being precisely the shot-noise term. 

The second term, quantified by /„ from Eq. (fT7j) . 
marks the excess noise when the radiation field is not 
in a coherent state. As mentioned above, /„ = for co- 
herent radiation, and f n = 1 for thermal radiation. Even 
when /„ is small, any nonzero /„ will eventually become 
the dominating effect for the photonoise as its contri- 
bution increases quadratically with time, in contrast to 
all other contributions. The cross-over point, where this 
term becomes larger than the shot-noise term, is easily 
estimated from Eq. (fl8)) as r)Tf n I n > 1. For any ther- 
mal radiation this regime thus is already entered once 
the measurement time is long enough for one photon per 
mode to be counted on average. I n can be computed from 
the Bose-Einstein factor, evaluated at the temperature 9 
of the source, 



b(u,e) 



i 



exp(Htjj / 1cb8) — 1 



(19) 



This value depends on the ratio of frequency and tem- 
perature, and thus has the same value b sa 0.004 for both 
a light bulb emitting visible light and an EUV plasma 
source operating at up to 200 000 K to emit radiation 
with A « 13 nm. 

The third and final contribution in Eq. (|18p quantifies 
the beating of different frequency components when the 
measurement time is finite. It is sometimes referred to as 
dynamic laser speckle since it is most relevant for mea- 
surement times T not much larger than the coherence 
time r of the radiation, 



dt\G(t)\ 2 = / 6u\G(u)\ 



(20) 



hence the name "temporal degree of coherence" fre- 
quently given to G(t). However, the factor [T— \t\] is also 
relevant, and it would be wrong to replace the integral 
simply by J exp(— i/r) 2 dt, as we will now demonstrate. 
The extra factor [T — \t\] , which seems to be largely ig- 
nored in phenomenological literature, gives less weight to 
the frequently surprisingly large tails of G(t). 

Two frequently encountered spectral distributions are 
Gaussian and Lorentzian, 



G 



Gauss 



G 



Lorentz 



(u) = %/2t 



-2ttt z 



2t 



1 + 47T 2 T 2 W 2 



(21) 
(22) 



The curves above are already normalised to 1, i. e., their 
integral (not their square integral) is 1. The widths of 
the spectral distributions have been chosen such that the 
coherence time, computed from Eq. (|20[) , is the same 
value r, allowing for a direct comparison of these two 
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FIG. 1: Value of the integral over time in Eq. (|18|l that quan- 
tifies the strength of the dynamic speckle. The two curves at 
the bottom follow by inserting the expressions for a Gaussian 
and a Lorentzian spectrum from Eq. (I24|l . When the term 
(T — \t\) in Eq. (|18[1 is replaced by T, and an exponentially 
decreasing G(t) is assumed, the third curve follows. 



curves. Their Fourier transforms are 



GGauss(t) 
GLorentz (t) 



1 



/27T 
1 

7^ 



exp 
exp 



8ttt 2 

J! 

2ttt 



(23) 
(24) 



The results for entering these expressions into the integral 
in Eq. (|18[) are shown in Fig. [TJ There is a significant 
dependence on the spectral shape, with the longer tails of 
the Lorentzian leading to a slower approach to saturation. 
Please note that the graph corresponds to the variance of 
the photocount. When the square root of the variance is 
scaled by the mean intensity - this quantity is called the 
contrast - the saturating curve shown in the figure turns 
into a curve decreasing to zero as T becomes larger. 

Gxorentz is a plain exponential and thus conforms best 
to the simple model of how to model finite temporal co- 
herence. We thus have used this curve for contrasting 
Eq. (|18|) by an approximation where the factor T — \t\ 
is replaced by T, thereby ignoring the additional depen- 
dence on \t\. For all values of T, the integral increases 
by about 40%, making the factor T — \t\ important for 
correct modelling of the photon counting statistics. 

A final word about the cross-over from "dynamics" to 
"statics" seems to be in order. Since the temporal effect 
begins to saturate at time r it is frequently assumed that 
the nontemporal regime is entered once T 3> t. While 
this statement is correct for single-mode radiation fields, 
it is incorrect here since the shot-noise term n, setting 
the reference value, scales linearly with the number N of 
modes whereas the prefactor of the temporal term scales 
as N ■ (N — 1) such that the condition T ^> t has to be 
replaced by T ^ Nr. Incidentally, since it is the regime 
Nt »T>t that current high-power excimer lasers are 
operating in Q, this distinction is of actual importance. 
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IV. ETENDUE 

Similar to classical mechanics, theoretical optics knows 
the concept of phase space [Io| • The optical phase space 
is spanned by position and spatial frequency, informally 
called fc-vector, and the radiation field is completely de- 
scribed by the (pseudo)-density function W(f, k) known 
as the Wigner function. Any radiation field then occupies 
a certain volume in phase space. Unfortunately, there is 
no good mathematical metric to actually measure or at 
least define such a volume. 

In experimental and technological areas of optics, in 
contrast, the volume of phase space occupied is well- 
defined and is referred to as etendue. The energy density 
I(f, a) of the light, given as a function of position and 
angle, is assumed to be larger than zero only inside a 
finite and well-defined region. This allows for an easy 
and direct definition of the etendue of the radiation field. 
Etendue is an important quantity since, at least within 
geometrical optics, it can only be changed by vignetting 
the radiation field. The etendue supplied by some radia- 
tion field incident on some optical apparatus should thus 
be no larger than the etendue that this apparatus can 
accept. 

Convolving the Wigner function W(f, k) with a small 
kernel, its size given by an uncertainty-relation [llj . 
yields I(r, a) when the relation 



a = kX , 



(25) 



between angle and fc-vector is utilised. Thus on first sight 
it might seem that the relation between microscopic and 
macroscopic definition would be straight-forward. 

The essential difference is that W(r, k) has no finite 
support whereas I(r, a) is assumed to have. It can be 
shown that due to the uncertainty principle, every mode 
<fi n needs to have infinite tails in at least either real space 
(position) or fc-space (angle) [12j . and this is reflected in 
W(f, k). Of course, also I(r, a) has the same infinite tails 
but in the macroscopic world, the function values inside 
those tails drop down so quickly that they are "assumed" 
to be zero once they are so low that they have become 
irrelevant for practical purposes. 

We thus need a definition for phase space volume that 
is based on a microscopic radiation field but reduces to 
the macroscopic etendue for "large" radiation fields. For 
the purpose of this paper, the relevant way of quantifying 
phase space volume is by the number of modes that can 
be fit into it. We will tackle the problem in the opposite 
way: We pick a certain number of "sensible" base func- 
tions and compute all possible radiation fields that can 
be formed from them. We will then identify for the in- 
tensity profiles I(r) as function of position and 1(a) as a 
function of direction the properties that allow to deduce 
the number of base functions used. 

To have shorter mathematical expressions, we will ex- 
press directions in terms of the fc-vector instead of the 
angle a. Furthermore we make two assumptions. First, 



we restrict ourselves to one spatial dimension but the ex- 
tension is straight-forward and becomes trivial when the 
radiation field factorises in its spatial dimensions. Sec- 
ond, we assume that the intensity distribution I(x, k) 
factorises into a spatial term I^ x > and an angular term 
l( k \ Since the etendue accepted by practically any real- 
life apparatus factorises, these assumptions do not limit 
the applicability of this approach. 

We start with N + l orthonormal base functions <fi n {x), 
n = 0, . . . , N. The orthogonality is essential since this 
prevents the base functions from being too similar and 
thus collapsing into a very small volume of phase space. 
Furthermore, all modes of the radiation field have to be 
mutually orthogonal anyhow as they are eigenfunctions 
of a Hermitian operator, cf. Eq. @. We assume them to 
be normalised for the ease of the calculation. For every 
(fi n (x), there also exists its Fourier transform 



1 



<fi n (x)e ikx dx 



(26) 



To progress, we need to choose a certain set of base 
functions. The aim is to pack the base functions as 
tightly as possible in phase space, and without loss of 
generality we try to pack them around the origin. If all 
<fi n , n — 0, . . . , N, are assumed to fulfil the conditions 



\<fin(x)\ <C(1 + \X\) 

\4> n (k)\<C(i + \k\) N 



N c ~a\x[ 



-b\k\ J 



(27a) 
(27b) 



then it can be shown [l3[ that ab = 1/4 is the most 
condensed situation for which there are solutions. Fur- 
thermore, all solutions can then be written as 



(fin(x) = P N (x)e- a ^ 2 



$ n (k) = Piv(fc)e- b|fe|2 (28) 



where Pn(x), -P/v(fc) are polynomials of order not higher 
than N. Only the product ab is fixed by this but not a or 
b on their own. This freedom is equivalent to changing 
(finix) — > (fi n {x/r) with the simultaneous change <fi n (k) — > 
(fi n (kr). This also reflects that, while phase space volume 
is well-defined and preserved as volume, its projection 
onto a particular axis is allowed to change. 

The set of solutions (fi n (x) of Eq. (f28|) is N + 1 dimen- 
sional, and a convenient orthonormal set consists of the 
Hermite functions ifi n (x), defined as 



ifin(x/r) 



1 



■. exp 



2r 2 



H n {x/r) , (29) 



where H n (x) are the Hermite polynomials, 

H (x) = l, H 1 (x)=2x, 

H n+ i(x) = 2xH n {x) - 2nH n -x{x) . (30) 

The prefactor in Eq. (|29|) ensures normalisation of ifi n (x). 
The Hermite functions are eigenfunctions of the Fourier 
transform, 



(31) 
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such that the derivation presented in the following for 
real space also applies to fc-space. 

All normalised functions (j>(x) inside the solution space 
of Eq. (f2"T)) can be written as 



N 



(X) = ^ a nlpn(3 



N 



= 1 



(32) 



The maximum intensity |(?!>(a;)| 2 possible for such a state 
can be computed using the method of Lagrange multipli- 
ers. Finding an extremum of <fr(x) is equivalent to finding 
one of |<Tj(a;)| 2 , and we will pick either of these quantities 
depending on which one is more convenient. The con- 
dition dA/dao = . . . = dk/daN — dA/dfi = of the 
function 



JY 



N 



A(oo, ■ ■ ■ , ajv, m) = X! a n' l l ) n{x) - - ^|a„| 2 ) , (33) 



n=0 



n=0 



yields the intermediary result ipk(x) + 2/ia.k = 0. By 
multiplying this with tpk(x) respectively at and summing, 
this gives the two conditions 



JY 



N 



J2Hn(x)\ 2 = -2n^a n tl) n (x) = -2n<j>(x) , (34a) 

N N 

-2^\a n \ 2 = ^a n ^„(a:) = <j>{x) . (34b) 



n=0 



n=0 



Solving for fi by eliminating 4>(x) and inserting /z into 

Eq. (|34bj) . remembering that X)n=ol a "l 2 = •"■> one arr i ves 
at the result for the maximum intensity, namely 



N 

rnaxl^)! 2 = Y,\Mx)\ 1 



(35) 



The method of Lagrange multipliers only gives necessary 
but not sufficient conditions for extrema, meaning that 
there cannot be any additional extrema different from 
Eq. (|35| but this equation might not describe an ex- 
tremum in the first place. Equation (|35[) actually de- 
scribes two solutions, <j)(x) = + . . . and <j){x) =—.... It 
is obvious from physics reasons that at least two extrema, 
one being the most positive and one the most negative, 
exist. Hence, Eq. ([33)) uniquely describes these two ex- 
trema. 

A graphical display of the solutions max|</>(x)| and 
max|<j!>(a:) | can be found in Fig. [2] as a function of TV. The 
right edge of max|0(x)| becomes progressively steeper as 
TV increases. This means that in the macroscopic limit 
TV ^> I there exist a sharp edge that can be associ- 
ated with the macroscopic etendue. To define the po- 
sition of this edge, one could pick the ^-coordinate where 
max|</)(x)| 2 = 1/2 or where it is equal some other particu- 
lar value but this would lead to an arbitrary result. It can 
be shown that the number of basis functions is then pro- 
portional to the square of the width at which the "cut" 




FIG. 2: Maximum values max 1 4>(x) and max|<^>(fc)|, computed 
from Eq. (|35[) , as a function of TV where TV is the number of 
basis functions minus 1. Depicted are the curves for TV = 
1, 4, 9, ... , 13 2 . If the corresponding curves for a normalised 
radiation field are below a given line, the radiation field cannot 
have more modes than indicated by the number TV on the 
curve. 



is done [14| such that the result depends more on this 
arbitrary decision than on the radiation field. Also the 
method of "almost bandwidth-limited functions" (which 
is usually formulated for time-dependent electric signals, 
hence the term "bandwidth" instead of "angular range" ) 
suffers from a similar problem [l5j |. 

Luckily, there exists a strong link between the micro- 
scopic and macroscopic worlds because the Hermite func- 
tions tpn(x) are the solutions of the eigenvalue equation 



d 2 ip n (x) 
da; 2 



+ (2n+l-a: 2 )V>„(a;) =0 



(36) 



that can be found in any quantum mechanics textbook 
as it describes the harmonic oscillator. The eigenfunction 
ip n has eigenvalue (=cncrgy) of n + 1/2, and a classical 
oscillator with the same energy is restricted to the in- 
terval —y/2n + 1 < x < \j2n + 1. The extremal curve 
4>(x) computed above thus is the extremal superposition 
of harmonic-oscillator solutions with energy level n up to 
N. Hence, the classical edge beyond which the oscillator 
cannot be found is given by y/2N + 1. 

Summarising, a microscopic radiation field has at most 
N + 1 modes when for all of its modes <j> m (x) the condi- 
tions 



jV 

\^rn{x)\ 2 < Y,\MxMf 
n=0 
N 

l<Ufc)l 2 < ]>>»( fcr )l 2 



(37a) 
(37b) 



are fulfilled for some value r identical for all modes. Since 
|0m(aO| 2 an< i I0m(fc)| 2 are the intensities as a function of 
position and angle, respectively, these are quantities that 
can in principle be measured. 

A macroscopic radiation field with a half-diameter X 
in real space and a half-diameter K in angle space is thus 
spanned by modes ip n with n = 0, . . . , N determined by 



V 2TV + 1 = X 



1 



V2N + 1 = K 



(38) 
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The scale factor r drops out when expressing results in 
terms of the macroscopic ctenduc. For given etendue 
E = XK, the number N of modes (minus 1) is thus 



S' 



N 



1 



(39) 



If the etendue is specified in angle (radians) instead of a 
k- vector, this becomes 



N 



E/X - 1 



(40) 



STATIC SPECKLE 



The variance of the photocount computed in Sec. IIIII 
describes the uncertainty of the photocount if the ra- 
diation field at the detector is completely known. Fre- 
quently, the mode structure is known at the light source 
but not at the detector. The standard undergraduate 
example is a HeNe-laser pointed at the wall of the lec- 
ture room where the observer can see a speckle pattern. 
When the eyes are then moved, the direction of motion 
of the observed speckle pattern depends on whether the 
viewer is near-sighted or far-sighted The HeNe-laser 
emits a very regular mode structure but the modes are 
then changed by scattering at the rough wall. The influ- 
ence of the eye movement demonstrates that scattering 
at the wall does not result in local intensity changes at 
the place of scattering but rather in a deformation of 
the mode structure that translates into intensity changes 
only by propagation to and focusing in the eye. 

The change of the mode structure on its way from the 
radiation source to the detector can be unintentional, just 
as in the example above, or it can be intentional due to 
some apparatus designed to mix the modes of the radi- 
ation field. Unless optical systems are manufactured to 
the highest standards, they will have surface roughness 
in excess of the wavelengths. This results in uncontrolled 
- thus in a certain sense random - mode deformations. 
In addition, any apparatus designed to provide a uni- 
form and stable illumination at its exit can only do this 
by means of "mixing" the incoming light. Illumination 
systems used in optical lithography are the supreme ex- 
ample of such an apparatus as the local intensity at the 
output varies only by about 1 % even when the form of 
the illumination of its entrance is completely changed. 
While the light paths inside such an apparatus might be 
very controlled, for an outside observer there (intention- 
ally) is no recognisable connection between the light at 
the entrance and the exit. In other words, an apparently 
random change of the mode structure occurs in this case 
as well. 

In Sec. IIIII we have used the modes </> ra (^) of the radi- 
ation field to expand the operator for the electric field, 
cf. Eq. (|4]) . Any orthonormal basis could have been used 
instead but using the modes offered the advantage that 
the associated annihilation operators a n {uj) become un- 
corrected, cf. Eq. 0. We will now switch to a different 



basis transformation 



light source 



«=i,...^v 



S" 



"apparatus" 
to deform 
mode structure 



FIG. 3: Graphical display of the transformation steps to ar- 
rive at static speckle. By means of a basis transformation, 
described by S' , one switches from the modes <f> n that are 
optimal to describe the light source, to new base functions 
<f)' n that are optimal to describe the photodetector. Along 
with the transformation tj> — > </>', the annihilation operators 
also need to be transformed, a —> a' . The additional action 
of scattering or mixing is then described by the matrix S" , 
transforming a' into 6. 



basis set, at the expense of replacing Eq. ([7]) by some- 
thing more complicated. 

Considering only the radiation source and the detec- 
tor, hence ignoring the scattering between them for the 
moment, we pick new basis functions <j)' n {r) and associ- 
ated annihilation operators a' n such that at the location 
R of the detector all 4>' n except for <p[ vanish. Such a 
basis set can always be found and offers the advantage 
that photodetection at the point R has been reduced to 
the photocount of the operator a\. Since both bases are 
orthonormal, the transformation between a = a±, <Z2, . . . 
and a' = a'±, a' 2 , . . . can be described by a unitary matrix 
S', 



i = s'c 



(41) 



with SS^ = 1 . Computing S' explicitly would be a very 
formidable task but it will turn out that knowledge of >S" 
is not needed. 

We will fix the number M of new basis functions <fi' n 
shortly but need to allow for the case that M is larger 
than the number N of excited modes of the radiation 
field. This is easily achieved by simply adding M — N 
vacuum states to the input for Eq. (141]) . 

We now assume that the radiation field is perturbed 
between the light source and the place of the detector by 
means of some "virtual device" - either intentionally by 
a properly designed apparatus or, usually unintention- 
ally, by quasi-random scattering. This can be described 
in two opposite but physically equivalent approaches. Ei- 
ther, one assumes that the mode structure is perturbed, 
thus changing {<fi' n } while keeping {a' n } unchanged, or 
one assumes that the mode structure {4>' n } is unchanged 
while energy is transferred between modes, thus modify- 
ing {a' n }. We use the second approach, also known as the 
method of input-output relations , since the annihila- 
tion operators {b n } describing the radiation leaving the 
device can then be expressed in terms of the annihilation 
operators {a^} entering device. In the absence of non- 
linear media, this relation is linear, and can be described 



8 



by a matrix S" 



b = S" 



(42) 



If there is no absorption, S" is unitary. Apart from en- 
ergy conservation, already the commutation relation Q 
demands the unitarity of S. In the presence of absorp- 
tion, Eq. fl4"2")) would need to be supplemented by noise 
sources [16( to ensure these commutation relations. 

In the absence of further information, the unitary ma- 
trix S' describing the distortion of the mode structure 
is uniformly distributed in the space of unitary matri- 
ces [l7j | . This concept of "maximum uncertainty" is 
known from other areas of physics as well, most notably 
from thermodynamics where the basic lemma is that ev- 
ery microstate compatible with the macroscopic bound- 
ary conditions is equally likely. 

The photocount distribution at some point R has thus 
been reduced to the photocount distribution for the mode 
b\ , related to the original incident mode structure of the 
radiation source by 



S3 , 



S = S'S" 



(43) 



Since for a given position of detector and light source R, 
S' is a fixed unitary matrix, the product S = S'S" is also 
distributed uniformly [l8|], and so no knowledge of S' is 
needed, as promised above. 

The only information about what happens between 
light source and detector enters in the form of the number 
M of basis functions. This number is not arbitrary but 
is determined by the physics of the mode deformation 
process described by the matrix S". Please remember 
that M is the number of base functions that are mixed 
by the transformation S" . If the mode deformation is 
due to some intentional mixing in an apparatus accept- 
ing and emitting some etendue E, the number M follows 
from the results of Sec. IIVI and has a finite, well-defined 
value. If the mode deformation is due to scattering or a 
diffuscr, M is determined by the number of basis func- 
tions that could in principle end up in the detector due to 
scattering or diffusing. Depending on maximum scatter- 
ing angle, this can be a very large number and is always 
larger than N. 

Computing the first factorial moment analogue to 
Eq. JIT]) yields 

n« = n( [ J dtb\(t)h(t)) 



nS nl )I n 



is equal to 1/M where M is the order of the matrix, 
yielding 



where in the final equality sign all averages over the fluc- 
tuations of electromagnetic field have been taken, just as 
described in Sec. lIIIl and an average over the ensemble of 
random unitary matrices remains. The average is easily 
computed from SS^ = 1, hence the average in Eq. (14"4"1) 



M 



(45) 



This is the expected result, namely that the mean inten- 
sity at the detector is equal to the properly scaled total 
intensity of the original light field. 

Computing the first moment demonstrated that to 
translate the equations from Sec. IIII1 all terms of the 
form 4> n (R) have to be replaced by Si n , and all terms 
4>n{R) by S} lX . The correct starting point is Eq. (fl"6j) . 
which transforms into 



,(2) 



V 2 T 2 



da;^(|5 ln | 4 )([at( W )] 2 [a n (a;)] i 



V 2 T 2 £ (\SmA 2 \Sin 3 \ 2 )I n Jn 
„2 r T 



J dt £ (|S lni | 2 |S lT 



rj_ 

2ir j _t 



x I ni I n2 G ni (t)G* n2 (t)[T - \t\] . (46) 

It only remains to compute the necessary averages over 
the unitary group in a similar spirit to above where we 
found that (SinS^) = 1/M. The average of the square 
of this expression is frequently needed and can be found 
in many papers, e. g. in Ref. 19, with the result 



{[Si n S nl ] ) 



M(M + 1) 



(47) 



The average (|«S'i n | 2 |S'i m | 2 ) for n 7^ m is not identical 
to (|<S'in| 2 )(|5'i m | 2 ) as the unitary condition SS' = 1 in- 
troduces correlations among the elements of S. Starting 
from the exact relation 1 = S ra ml'^ lm l 2 |^' llI l 2 auows us 
to write 

E(i^i 2 i^i 2 )+E(i^i 4 ) = 1 ' ( 48 ) 

which, together with Eq. (|47p . gives us the desired aver- 
age, 



M(A/+1) n ¥ zm 

M(A/+i) n = m 



(|Si„| 2 |Si 

This allows us to transform Eq. (|46j) into 

n 2 T 2 
12) V * 



(49) 



(44) 



2 ,; (tf + 1) j^E([4Hl 2 KMf) 

+ M(M + 1) 4- ni " 



2ttM(M + 1) 



T)f_ T dt E InJuAnWUtW-W 



(50) 
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When computing the variance in Sec. IIII1 the term with 
Iniln 3 dropped out when [ro^] 2 was subtracted from 
n^ 2 \ This is no longer the case here, and one arrives 
at 



1 



var n = n - 



-n 2 



ifT 2 



(M + l) M(M+1) 



(51) 



Before discussing this result, we will analyse an approx- 
imation that is frequently done: The exact averages from 
Eq. (I49p are replaced by their Gaussian approximation, 
equivalent to a phasor representation 



|S'ln| 2 |S , lm| 2 ) 




n 7^ m 
n = m 



(52) 



The denominator in Eq. |g5J) was M(M + 1) so that 
the mixing of the modes would conserve energy. Equa- 
tion (j52")l neglects this constraint. While this error might 
seem to be negligible for large M, we will demonstrate 
that its effect is surprisingly large. First, we compute 



,(2> 



rfT 2 
2n 2 M 2 



rfT 2 
If 2 



V / i 

/ i ni i n 



V 



2ttM 2 



pJ_ T dt E I n Jn 2 G ni (t)G* n2 (t)l T -\t\] 



(53) 



In contrast to Eq. (|50j) . the term with I ni In 2 now cancels 
against [tt/ 1 )] 2 . The variance then becomes 



„ 2 T 2 

vav».»li+-^23(/ n + l)J* 



M 2 



>1 



2ttM 2 



(54) 



This is basically Eq. (jT5J) with an additional contribution 

V 2 T 2 



var n sta tic 



M 2 



£3 



(55) 



that is called "static" as its scaling with T 2 implies that 
its strength, when scaled by the mean squared intensity, 
is independent of measurement time T. 

We arrived at the ensemble of random S for calculating 
the variance var n sta ti c by assuming a fixed position of the 
detector while the mode deformations between radiation 



source and detector is different for a every member of 
the ensemble. However, from Eq. (|4"3"|) it follows that 
this variance can equally well be computed or measured 
by keeping the mode formations constant while changing 
the position of the detector. The static speckle contrast 



c, 



= yvar n sta tic _ [En /2 J 1/2 



(56) 



can thus equally well be determined by doing a (long) 
measurement at different detector positions R. In prac- 
tise, one uses a camera instead of a point-like photode- 
tector, and c is calculated from the observed contrast of 
the image taken by the camera. The image resembles the 
speckle pattern found on the fur or skin of many animals, 
hence the name "speckle". The static speckle contrast 
can, in this approximation, conveniently be expressed in 
terms of the effective number of modes N c g, 



N eS = 



[En Inf 

y i 2 



(57) 



that is equal to the actual number of modes if all intensi- 
ties are equal, and smaller (implying more static speckle) 
otherwise. 

For the exact result ([5T|) , static speckle is not so easy to 
define, as the difference between Eq. (f5"Tj) and the results 
from Sec. IIIII is not a simple term proportional to T 2 . 
However, taking the limit T — > oo, hence living up to the 
name "static", the difference between these two terms 
gives 



[ var "Eq.tSB 



var n 



Eq 



1/2 



M 



f 



M+1 [EM 

M 



M 



-fJ- 

1 ViVeff 



1 

M 



M+1 [E n i n y 

, M-iy n fnll 



M + l 

1/2 



1/2 



M + l 



[En In] 



(58) 



There thus is a very significant difference between the 
exact result and the approximation. In the latter, the 
contribution (|55p is minimal but still nonzero if all modes 
carry the same intensity /„ whereas for the exact result 
Eq. (|5ip . assuming coherent radiation, this contribution 
vanishes if all /„ are identical since then N e g = M, 

The approximation becomes valid in the limit that 
there are many modes but only a small number of them 
are actually excited: Equation (|55|) depends on the av- 
erage squared intensity (I 2 ) n of the modes of the light 
source, where for a fixed light source the average is taken 
over its modes. This kind of average is indicated by the 
notation (...)„. In contrast, the exact result Eq. (|5l"j) 
depends on M/(M + 1)[(/ 2 )„ - (I„) 2 ] but the term (J n ) 2 
scales with the square of the fraction of excited modes 
and may thus be neglected in the limit stated above. 
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This important conclusion can also be formulated in 
terms of the etendue introduced in Sec. IIV1 The approx- 
imation is valid in the limit that the light source fills only 
a small fraction of the etendue that the "device" used for 
mixing can accept. It also means that mixing will not 
introduce static speckle as long as the entire etendue of 
the mixer is uniformly filled by the light source (assuming 
coherent radiation). 

This might seem to contradict phenomenological theo- 
ries that claim that there should aways be some pattern 
in the detector plane because light from the same point 
of the light source can reach the detector along differ- 
ent paths, thereby resulting in uncontrolled interference 
that inevitably creates intensity fluctuations. However, 
the unitary condition SS* = 1 ensures that, whenever 
a mode of the light source is brighter at the detector, 
another mode has to be dimmer. This property is ful- 
filled for any individual S since Enl^i"! 2 = ^ an( ^ 110 
averaging over S is necessary to arrive at this conclusion. 

Technologically, mixing apparatuses have a well- 
defined etendue, hence a well-defined value of M. When 
the etendue that they can accept, is filled to a large ex- 
tend by the incident radiation, the effect of the unitary 
condition is strong, and the exact formula (|5l"j) has to be 
used. For random scattering, on the other hand, in par- 
ticular by a diffuser, frequently N <C M, and therefore 
the approximate formula (|55[) might suffice. 



VI. STATIC SPECKLE AND ETENDUE 

In Sec. |V] we demonstrated that the static speckle de- 
pends on the mean squared intensity in each mode of the 
radiation source, hence on J^n ^n> c f- Eqs. (fB"T) and ([53)) . 
The effect is minimal if all modes carry equal weight, 
and the maximum number of modes possible follows from 
the etendue via Sec. II VI The minimum amount of static- 
speckle possible can then be computed - it is zero if the 
etendue is not increased by the scattering and mixing 
processes, and finite otherwise. 

This raises the following question: is it possible not 
only to compute the minimum amount of static speckle 
for given etendue of the radiation field but also the "typ- 
ical" amount of speckle? We will demonstrate that static 
speckle contrast is a self-averaging quantity, i.e., in the 
limit of many modes M, the relative fluctuation while 
taking the average over the unitary group becomes small. 
This implies that the ensemble average is characteristic 
for almost all realisation of the ensemble. 

When cf> n (r) are the modes of the electric field, the 
electric field E(r) is uniquely described by coefficients 



(59) 



and the static speckle contrast is determined by the effec- 
tive number N c g of modes in a more or less complicated 



way depending on whether the exact result or the ap- 
proximation is used, 



1 



(60) 



In Sec. lIVI thc linear span corresponding to some given 
etendue E was computed. In the absence of dynamic 
information, one cannot select the modes out of all pos- 
sible sets of basis vectors, and any set of basis vectors 
could thus represent the modes. We assume that the de- 
composition (|5"9"|) represents the actual modes (f> n (x) of 
the radiation field. Since, as just explained, we cannot 
know (j> n (x), the equations for the static speckle contrast 
are instead evaluated using some other orthonormal basis 
<p' n (r). The two bases are related by a unitary transfor- 
mation, <p{r) — U<f>'(f), and consequently a? = Ua, but 
U is unknown for the reasons just explained. 

If no additional information is available, the best avail- 
able estimate is the average of 



1 



E 



n K\ 4 



N *e [E„Klf 

[E„M 2 ] 2 



(61) 



after integrating over the unitary group, i.e., over all 
possible unitary matrices. We study the average oil/N^ s 
instead of N^ s since it both is easier to compute and is the 
relevant quantity for the speckle contrast. The variance 
([1/A^ ff ] 2 ) - (l/iV^ff) 2 then specifies how well the average 
is representative for all possible members of the ensemble. 

Averaging over the unitary group, the nonvanishing 
terms are of the form 

( U il,h ' ' ' U iM,3M U P(il).P'(jl) ■ ■ ■ U P(iM),P'(j M )) 

= V s(P - lpi) , (62) 

where P and P' are permutations of the numbers 
1, . . . , M, and s denotes the cycle structure of a per- 
mutation. The coefficients V depend only on the cycle 
structure of P~ 1 P' and have been tabulated [l9j . 
Labelling the numerator in Eq. (|6"Tj) as Z, 



n,ti ,...,*4 



(63) 

and counting all the permutations of the indices in it 
yields 

(Z)=2(^, 1 + V 2 )^^| ail | 2 K| 2 . (64) 



For the square, 



^ 2 = EKI 4 K 



4 

m I i 



(65) 
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one finds 

(Z 2 ) = {W 1AAA +24y 2 ,i,! + 32^,! + W 2i2 + 24F 4 )x 

^ ^ ^ ^ 1^21 a «2 Q U I 

n^m zl,...,?4 

+ (241/i,i 4 ,i + 144^2,1,1 + 192F 3 ,i + 72V 2 ,2 + 144V 4 )x 

Y \ a Ha l2 a i3 a i4 \ 2 . (66) 

n zl,...,i 4 

The terms , . . . , Oi 4 appearing in the averages above 



cancel against the denominator of Eq. (fBTj) . This is as 
expected since, as explained, these coefficients depend 
on the base that was used to determine them, i. e., they 
depend on the choice of one particular unitary matrix. 
Averaging over the unitary group, this choice must be 
irrelevant, thus the coefficients drop out, and the results 
depends only on the number M of modes. 

Inserting Eq. into Eq. (fBTj) and using the tabulated 
coefficients from Ref. [lj^ gives the result 



M + l 



(67) 



which has an easy interpretation. The term M+l instead 
of the naive term M in the denominator is due to the 
correlations imposed on U as it is unitary, cf. the factor 
M + 1 in Eq. ([55)1 . The 2 in the numerator implies that 
every mode is occupied with an effective intensity of 1 /2 
relative to the maximum intensity of any mode - a very 
plausible result since the intensity is a random quantity 
between and that maximum value. 

A similar calculation using Eq. (|66|) gives for the square 

1 M + 5 (68) 



(Af+l)(M + 2)(M + 3) 

so that the variance becomes 

1 M - 1 



Var N^ s 4 (M + l) 2 (M + 2)(M + 3) ' ^ 

The variance decreases much faster with M than 1 /N' cji 
or l/iVg ff 2 as the lowest-orders in 1/M cancel when com- 
puting the variance. This implies that l/N' eS is a self- 
averaging quantity, and its value computed from the 
number M and hence from the etendue of the radiation 
field gives a good estimation of the actual value of l/N' eS . 

One should be aware, however, that this is a statistical 
statement: in "almost all cases" this statement is correct. 
By use of particular light sources it still is possible to have 
a very different value for l/N' cS . This is also reflected in 
the calculation presented here. Equation (|6"0|) gives the 
actual value of 1/N e g which, of course, depends on the 
values of the a n . The average value from Eq. (f67|) no 
longer depends on the a n and can thus deviate. 

VII. NONORTHOGONAL MODES 

The modes <p n (r) of the electromagnetic field are mu- 
tually orthogonal and stationary, cf. Eq. ©, and all the 



time-dependence of the electromagnetic field is included 
in the frequency-dependence of the annihilation opera- 
tors. However, frequently it is rather inconvenient to de- 
scribe time-dependence in this way. An extreme example 
are pulsed lasers where it is more appropriate to describe 
the mode structure of each pulse separately. This comes 
at a price, though, as the modes of one pulse are then 
not necessarily orthogonal to the modes of another pulse. 
This is primarily an issue for the static speckle computed 
in Sec. [V] as the mutual orthogonality of the modes was 
essential for the derivation presented there. 

The static speckle contrast in Eq. ([58")) has an easy in- 
terpretation: except for a minor correction due to the 
unitary of the scattering matrix, every mode creates a 
mutually uncorrelated speckle pattern in the detector 
plane. The effect of a superposition of different modes 
is easily understood from this, and for the static part it 
is irrelevant if all the modes are excited simultaneously or 
rather one after the other. In this section, we will show 
that the speckle patterns of two nonorthogonal modes 
(i.e., necessarily at two different times) are correlated. 
We will, without loss of generality, refer to these two 
times as first and second pulse, respectively. 

The first pulse consists of modes 4> n and annihilation 
operators a n . For a fixed arrangement of mode structure 
of the light source, scattering between source and detec- 
tor, and detector position, thus only averaging over the 
fluctuations of the radiation field, the mean photocount 
is, completely analogue to Eq. (|4~X| . 



(70) 



We now consider a second set (j>' n of modes of the light 
source but keep the scattering between source and detec- 
tor and the detector position fixed. The transformation 
from the 4>n to the <j>' n is linear and can thus be described 
by a matrix S' , 



S" 



(71) 



where the notation (. . . | . . .) denotes the overlap integral. 
The annihilation operators transform as 



k ' 



and the mean photocount thus becomes 
n' = vTY(SS')in(SS')i 1 i: i 



(72) 



(73) 



The essential quantity is the cross-correlation between n 
and n' after averaging over all possible scattering config- 
urations, hence over all S 1 , and only the matrix S' may 
remain in the result. A positive correlation means that 
the light patterns caused by the two radiation fields are 
similar. When the photocounts of two such fields are 
subsequently integrated on the same detector, the static 
speckle contrast will be higher than naively expected by 
the inverse square-root of number-of-modes law. 
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The cross-correlation between n and nl is given by 

c:={nn') - (n)(n') , (74) 

and the average has to be taken over the matrix S keeping 
the matrix S' fixed. Writing out the first term explicitly, 



nn 



nrnkl 



SinSi n SikSuS km S lm I n I m , (75) 



the average yields nonzero terms only for k = I, and the 
two cases k — n and k ^ have to be distinguished. The 
necessary averages have already been given in Eq. (|49|) . 
and the result is 



{nn 



M(M + 1) 



El q' \2t jl 
\ u nm\ J n J m 



^EE^i 2 ^;. ( 76 ) 



M(M + 1) 



This can be simplified by noting that J2 n ^k\^'k 
1 — |S^ m | 2 . The cross-correlator then becomes 



' 12 



M(M 



I y 



\{<t>nWm)\ 2 - Tl 



M 



I I' 



(77) 



This equation gives the expected result in the two ex- 
treme situations. If all intensities are equal, the cross- 
correlation becomes zero since, as explained earlier in 
the context of Eq. (|58|) . the static speckle itself van- 
ishes, and thus also the cross-correlation does. If only a 
single mode n is excited, the cross-correlation becomes 
proportional to \{4>n\<t>' n )\ 2 ~ l/M. The average value 
of | (<^>jt, |0^t,} | 2 when assuming a random S' is equal to 
l/M, and the cross-correlator is again zero. If <fi n and 
4>' n are identical, the cross-correlator becomes maximal. 
The cross-correlator thus correctly quantifies if 4> n and <f>' n 
are stronger correlated than in a random configuration, 
and the overlap integral |(0n|</>m)| 2 is directly related to 
the static speckle pattern. 

When a photodetector sums over N pulses, the vari- 
ance var n of the integrated photocount follows from the 
variances varn^ of the fc-th pulse, given by Eq. (|5ip . and 
from the cross-correlator cjm, given by Eq. (|77p . between 
the fc-th and the l-th pulse, 



var n = E var n k + E Ckl 

k k^l 



(78) 



Apart from the rather obvious relation (|75|) . the re- 
sults from this section have another application. Fre- 
quently, radiation fields are expanded not in their modes 
but rather in a set of "convenient" functions, such as 
Gaussian beams or plane waves. Drawback of this ap- 
proach is that this set then usually has to be overcom- 
plete, hence it has more elements than there had been 
modes, and not all its elements are mutually orthogo- 
nal. For example, two Gaussian beams have an overlap 



that decreases exponentially with separation but is al- 
ways larger than zero. Two plane waves restricted to 
a finite interval have an overlap that decreases exponen- 
tially with the angle between the two waves but is always 
finite. 

Depending on the way that this expansion is done, 
some coherence information might be lost, namely when 
doing the expansion such that only the time-averaged in- 
tensities at all field-points are matched, and each <fi n is 
then assigned an intensity I n . This might sound like a 
bad approximation but effectively amounts to using the 
assumptions from Sec. IVII Inserting the /„ into Eq. (|5Tj) 
would yield a low static speckle contrast since the number 
of In's is large. However, the correct equation to use is 
Eq. (|78|) where every term is evaluated for a single func- 
tion (f> n , and the nonorthogonality of the <p n implies that 
the cross-correlation terms are large. This is as expected 
since introducing additional functions <f> n to describe the 
same system should not change its computed properties. 



VIII. DISCUSSION 

The purpose of this paper has been to give a con- 
sistent and concise description of multi-frequency and 
multi-mode effects on photodetection on a quantum level 
trying to avoid the use of uncontrolled approximations. 
Its focus lies on exactness, not on ease of presentation 
or application to some particular device. It thus comple- 
ments previous work, for example by Joseph W. Good- 
man or the group around Christer Rydberg [l|, Q . 

In Sec. IIIII we derived how the uncertainty of the pho- 
tocount depends on the measurement time and the tem- 
poral coherence function. Among others, it was shown 
that a description via the coherence time of the radiation 
is insufficient to arrive at an exact result. Similarly, in 
Sec. [V] not only previously-known approximations of the 
static speckle contrast were retrieved, namely an inverse 
square-root law with the effective number of modes, but 
in addition also effects of finite etendue accepted by the 
"mixing device" or of noncoherent, e.g. thermal, radi- 
ation followed from the chosen mathematical formalism 
without additional effort. 

Some of the results of this paper can be explained, al- 
beit only on a qualitative level, by the "brick wall" model: 
The radiation field is like a brick wall (horizontally: time, 
vertically: position), and only bricks in the same col- 
umn or row can cause interference effects. Counting the 
fraction of such bricks gives an intuitive explanation of 
the scaling of dynamic and static speckle contrast upon 
changing system parameters. If this is sufficient, the 
methods and results presented in this paper are a bit 
of an overkill. If, on the other hand, exact agreement 
between predication and performance of some device is 
essential, this paper offers the advantage of giving results 
using only well-defined input parameters, thus not rely- 
ing on effective parameters that need to be tuned until 
the desired result is retrieved. 
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Main function of any theory is to predict whether 
some device or experimental setup will work as expected 
and / or required. For state-of-the-art technological ap- 
plications, the performance might improve by about 25 % 
from one generation of a device or setup to the next - 
but rarely more. To predict if something actually is an 
improvement, the error due to approximations thus has 
to be much smaller this number. As an example, Fig. [1] 
gives an indication of the error made when using the co- 
herence time instead of the exact shape of the frequency 
spectrum. If improvement of dynamic speckle is a topic, 
one should thus use the exact formulae and refrain from 
approximations. 

Depending on the desired application, either the 
present paper or one of the previous works by other au- 
thors will thus be more appropriate. Furthermore, we 
restrict us to universal problems. For example, we study 
in Sec. fVl only the variance of the photocount, hence the 
contrast of the dark and bright "spots" in the plane of 
the detector, but not their spatial correlations, i.e., the 
size of these spots. Reason is that, while the contrast 
is universal, the size of the spots is not and depends on 
the setup. The book of Goodman [1] is to a large extend 
devoted to studying this question (on the classical level) 
for a large number of setups of technological or scientific 
importance. 

In Sec. lIIII the uncertainty of the photocount for known 
radiation field has been computed, yielding in addition 
to the well-known shot noise term two terms depending 
on the frequency spectrum and on the photon correla- 
tions of the radiation source, respectively. The impor- 
tance of the frequency spectrum is obvious, in particular 
in view of the two examples presented in this paper. For 
the importance of the effects of noncoherent radiation, 
the situation is less clear. For traditional thermal light 
sources, the necessary information can be found in this 
paper; for traditional lasers, the emitted radiation is co- 
herent, and the question becomes trivial. For modern 
light sources, such as excimer lasers |20| that are increas- 
ingly used for high-power applications, insufficient data is 
available. Technological advances are fast, and no good 
characterisation of the photon correlations for modern 
devices seems to have been published. Since this con- 
tribution to the noise inevitably becomes the dominant 



term once the measurement time is sufficiently long, this 
lack of knowledge is an actual issue. 

The static speckle, computed in Sec. E] depends on 
the number of excited modes of the radiation field (and 
the energy distribution among them) but, in addition, 
also on the number of modes that could in principle be 
excited. This is closely related to the concept of etendue, 
cf. Sec. lIVI Even though etendue is a very basic quantity, 
its microscopic definition does not seem to have been ad- 
dressed before. The relation presented between number 
of basis function and macroscopic etendue has applica- 
tion beyond this paper as any modelling of a radiation 
field - for the purpose of simulation or for an analytic 
study - needs to start with a given number of modes, and 
knowing this number a priori makes this process much 
more efficient. 

While the applicability of the relation between etendue 
and maximum number of modes is obvious, this is less the 
case for Sec. IVII discussing the relation between etendue 
and average effective number of modes. We demon- 
strated that this average is identical to the actual number 
for almost all radiation fields. However, that is "only" a 
statistical statement: it assumes that all possible radia- 
tion fields with given etendue are equally likely, and it 
is still allowed that a few cases (albeit of measure zero) 
deviate strongly. Different light sources have different 
properties, and for any technological application, one will 
choose the most appropriate light source. Selecting a 
single-mode laser that operates in a high-order Hermite- 
Gauss mode would give the most extreme disagreement 
between the actual number of modes (=one) and the pre- 
diction from Sec. IVII The results from this section can 
thus only be applied when the choice for a particular light 
source does not adversely impair the freedom of the mode 
structure of the generated light field. 
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